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Investigations of the polarization and relaxation of 209 Rn by spin exchange with laser optically 
pumped rubidium are reported. On the order of one million atoms per shot were collected in coated 
and uncoated glass cells. Gamma-ray anisotropies were measured as a signal of the alignment (second 
order moment of the polarization) resulting from the combination of polarization and quadrupole 
relaxation at the cell walls. The temperature dependence over the range 130° C to 220° C shows the 
anisotropies increasing with increasing temperature as the ratio of the spin exchange polarization 
rate to the wall relaxation rate increases faster than the rubidium polarization decreases. Polar- 
ization relaxation rates for coated and uncoated cells are presented. In addition, improved limits 
on the multipole mixing ratios of some of the main gamma-ray transitions have been extracted. 
These results are promising for electric dipole moment measurements of octupole-deformed 223 Rn 
and other isotopes, provided sufficient quantities of the rare isotopes can be produced. 

PACS numbers: 32.10.Dk, ll.30.-j, 23.20.En, 23.20.Gq 
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An electric dipole moment (EDM) is a T- violating sep- 
aration of charge along the spin of a system such as the 
neutron, an atom, or a molecule. Assuming CPT sym- 
metry, EDM measurements probe CP violation in the 
Standard Model and beyond. An atomic EDM would 
arise as a result of the presence of CP-odd fundamental 
interactions that affect physics at the atomic scale. Pos- 
sible examples include the vacuum expectation value of 
the gluon field (Oqcd), fundamental quark and electron 
EDMs, CP-odd four-quark interactions [lj], and the Wein- 
berg 3-gluon interaction [l], Q • The search for permanent 
EDMs has implications for many areas of beyond-the- 
Standard-Model physics. The expanded particle spec- 
trum introduced by supersymmetric theories increases 
the number of diagrams contributing CP- violating phases 
[|. The baryon asymmetry in the universe could be gen- 
erated by a mechanism that requires CP violation beyond 
that contained within the Standard Model [1] that could 
be revealed by EDM measurements. 

Isotopes featuring octupole deformation or vibrational 
strength, such as 223 Rn and 225 Ra, are expected to ex- 
hibit a large enhancement to their sensitivity to a CP- 
violating EDM [H, Calculations estimate that the 
octupole deformation of 223 Rn makes it approximately 
400 times more sensitive to CP violation than 199 Hg [5J], 
which currently sets the best limits on several mecha- 
nisms of CP violation. With a limit on the 199 Hg EDM of 
|d( 199 Hg)| < 2.10 x 10~ 28 e • cm Q, a measurement of the 
223 Rn EDM with a sensitivity of 10~ 26 e-cm would signif- 
icantly improve sensitivity to CP violation. NMR tech- 
niques have proven effective in making sensitive EDM 
measurements of some of the lighter noble gases [8[ , but 
they have never been applied to EDM measurements of 
small amounts of short-lived radon isotopes. In anticipa- 



tion of high 223 Rn production rates at TRIUMF's ISAC 
[9| and future facilities, a 209 Rn source was developed at 
the Stony Brook Francium Lab [l(| ■ Although the polar- 



ization of 209 Rn and 223 Rn has been demonstrated 111 ] . 



measurements had been made at only one cell temper- 

Rn 



ature. Both 209 Rn (1=5/2, T 1/2 =28.5 min) and 223 
(1=7/2, T 1/ / 2 =23.2 min), when polarized in a glass opti- 
cal pumping cell, are expected to have similar quadrupole 
interactions with the cell walls. Thus 209 Rn provides an 
appropriate system for the study of polarization and re- 
laxation processes. 

A beam of 16 incident on a gold target generated 
francium isotopes that were accelerated to 5 keV, elec- 
trostatically focussed, and implanted in a zirconium foil 
(see Fig. [TJ. Accelerating the 16 O beam to around 
91 MeV optimized the production, through the reaction 
197 Au( 16 0,4n), of 209 Fr (T 1/2 = 50.0 s), which has an 
11% branching ratio for electron capture decay to 209 Rn. 
The francium was implanted for about two half-lives of 
209 Rn. While there was some shielding between the foil 
and the detectors, the activity at the foil could still be 
monitored during implantation. For example, the 719 
keV line in the upper panel of Fig. [5] indicates 205 At in 
the foil, the product of 209 Fr a-decay (89% b.r.). After 
implantation, the radon was transferred to the optical 
pumping cell, which had been preloaded with rubidium, 
by simultaneously immersing the cell in liquid nitrogen 
and heating the foil to about 1000°C with the chamber 
isolated from the beamline and pumps. At this tempera- 
ture, radon rapidly diffused out of the foil [12]. Valve V2 
was then closed and about one atmosphere of nitrogen 
was added to the cell to act as a buffer gas, preventing 
radiation trapping. The cell was then isolated from the 
rest of the vacuum system by shutting valve V4, and a 
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FIG. 1: A schematic diagram of the experimental apparatus. 

two-piece glass oven was placed around it, allowing it to 
be heated to a desired temperature. Up to about one 
million 209 Rn atoms were transferred to the cell in each 
cycle. In the first of two runs, discussed in Ref. [HI 
and designated Run 1 in this letter, data were taken for 
an uncoated cell at temperatures ranging from 150 to 
200° C. In Run 2, data were obtained for an uncoated 
cell from 150 to 220°C and for a cell coated with octade- 
cyltrichlorosilane (OTS) Q from 130 to 180°C. 

The decay of polarized 209 Rn results in the emission of 
209 At gamma rays with an angular distribution depen- 
dent on the alignment, i.e. the second moment of the 
nuclear sublevel population distribution resulting from 
temperature-dependent polarization and relaxation pro- 
cesses. The cell was placed in a uniform magnetic field of 
about 10.5 gauss and illuminated by circularly polarized 
laser light from a Coherent diode laser system tuned to 
the Rb-Dl line. This polarized the rubidium, and spin- 
exchange collisions [151 ] transferred the valence electron 
polarization to the 209 Rn nuclei. The resulting angular 
distribution of gamma rays was monitored by HPGe de- 
tectors: a Eurisys Clover at 0° (4 crystals, each with 25% 
internal efficiency relative to Nal) and a 100% Ortec de- 
tector at 90° in Run 1, and two 100% Ortec detectors at 
both 0° and 90° (placed about two inches from the cell) 
in Run 2. The data from the individual crystals in the 
Clover detector were summed in a way that allowed it 
to be viewed as one large crystal. Data were also taken 
with the laser off in order to determine the unpolarized 



FIG. 2: Gamma-ray spectra from the 0° detector. Spectrum 
a) was taken over about an hour during the implantation of 
the 209 Fr beam in the Zr foil and b) was taken over the first 
ten minutes after transfer of the 209 Rn to the optical pumping 
cell. Note the different vertical scales. 

baseline rates at each detector and temperature. Earlier 
studies of the wall interactions of the spin-3/2 isotope 
131 Xe 16] indicated reduced quadrupole relaxation in un- 
coated cells, which motivated the study of both coated 
and uncoated cells in these runs. 

In order to extract information on polarization and 
relaxation, an optical pumping/spin-exchange model was 
developed to find the radon steady-state nuclear sublevel 
populations [171 ] . The rubidium polarization is calculated 
by numerical integration of laser-light absorption in the 
cell with the Rb-Rb spin destruction rate from Ref. [l8| 
and the rubidium number density given by [19J 

[Rb] = 10 9 - 318 -^Vfc B T, (1) 

following the method detailed in Ref. [2(|. At low tem- 
peratures, the rubidium is almost completely polarized as 
optical pumping overcomes all spin destruction processes. 
At higher temperatures, increased spin destruction with 
increased [Rb] leads to a decrease in the polarization due 
to limited laser power. The resulting radon polarization 
results from the interplay of processes whose effects can 
be parameterized by three main rates: the rubidium- 
radon spin-exchange rate 7s_e, the radon dipole relax- 
ation rate T\, and the radon quadrupole relaxation rate 
IV The spin-exchange rate is dependent on the spin- 
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exchange cross section use and the cell temperature T 
through the equation 



Ise = [Rb] (asEv) 



(2) 



in which the temperature dependence enters through 
[Rb] and the velocity-averaged rate constant (ctsev}. 
Quadrupole relaxation (r 2 ) has been shown to be much 
stronger than dipole relaxation (T\) in 131 Xe and 21 Ne 
H H3. Thus, ri is set to zero in the model. Wall re- 
laxation is modeled by an Arrhenius-type temperature 
dependence [22j], i.e. 



r 2 (T) = r?e To/T , 



(3) 



where fcsTo is the binding energy of the radon to the cell 
wall. The set of rate equations for the populations of 
the six nuclear sublevels of 209 Rn are written in terms of 
the rubidium polarization, jse, and r 2 , and an analytic 
equation for the steady-state sublevel populations is ob- 
tained. Selection rules are used to calculate the sublevel 
populations of the 209 At populated by electron capture. 
The angular distribution of the 209 At gamma rays de- 
pends on the moments of the populations, the spins of 
the states, and the multipolarity of the transitions [23] ]. 
Table Q] includes the multipole mixing ratios (6) of the 
four prominent lines and the spin transitions involved. 
Following the convention of Ref. [24|, 5 2 =a 2 /a|, where 
ai = 1 in a pure dipole transition, a 2 = 1 in a pure 




140 160 180 200 
Cell Temperature (°C) 



220 



FIG. 3: Anisotropy data from the coated and uncoated cells 
for the 337 keV 209 Rn gamma ray. For the uncoated cell the 
average of Runs 1 and 2 is displayed. R is the ratio defined 
in Eq. (4), and the solid curves are the fits, described in the 
text, from which a value for r§° is obtained. 
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FIG. 4: Anisotropy data for the 745 keV Rn gamma ray 
displayed in the same manner as the 337 keV data in Fig. [3] 
The solid curves are fits to the mixing ratio 8 given the value 
of r|° obtained from the 337 keV data in the corresponding 
cell type, as described in the text. 



1. The angular dis- 



quadrupole transition, and af+f 
tributions are expressed as a function of <jse, Tif, To, 
and 8 for each cell temperature under consideration. 

The 337, 408, and 745 keV lines were studied in de- 
tail. The 689 keV line could not be reliably extracted 
from the spectra due to the proximity of a gamma-ray 
line resulting from 209 At decay. The anisotropy for each 
gamma ray is given by the ratio 



R 



(n o /n 9 Qo) LON 
(uqo /ngo°) LOFF ' 



(4) 



where uqo /ngo° is the background-corrected ratio of 
counts in the photopeak at the 0° detector to counts 
in the 90° detector calculated for data with the laser 
on (LON) and with the laser off (LOFF). The 337 and 
745 keV lines show anisotropies that are consistent with 
those measured in Ref. The uncoated-cell 408 keV 

line (not shown here) shows no statistically significant 
anisotropy, also consistent with Ref. [111 ]. 

The solid lines in the figures were generated using the 
optical pumping/spin-exchange model. In Fig. [31 data 
for the 337 keV line are shown. Since <5 is known for this 
transition, the data can be used to extract the relaxation 
rate r£°, assuming values for use and To. For <jse, the 
estimate calculated by Walker [iH (2.5 x 10~ 5 A ) was 
used. An estimate for To (400) was made based on the 
131 Xe measurements reported in Ref. [lj|. For the 337 
keV data, best-fit values for rsf of 0.05 ± 0.01 Hz for the 



4 



uncoated cell and 0.032 ± 0.009 Hz for the coated cell 
were found. These fits are displayed as the solid lines in 
Fig. [31 There are multiple ways to generate temperature- 
dependent angular distributions that agree with the data 
using different values of <jse, To, and Tf, so the current 
data do not allow best fit values for all three parameters 
simultaneously. 

With the polarization and relaxation parameters fixed 
by the 337 keV data, the 745 keV and 408 keV data can 
be fit for S. The uncoated-cell 745 keV data allow (95% 
C.L.) values of 5 in the ranges <5 < 0.19 and <5 > 10.9. 
Previous measurements [26| exclude the range 6 < 2.86, 
and this work significantly improves that constraint. A 
similar analysis for the uncoated-cell 408 keV data allows 
the range 0.01 < S < 0.68 (95% C.L.). This establishes 
an upper limit of S < 0.68. The coated-cell data provide 
much weaker constraints. These results are not very sen- 
sitive to the values of <jse, To, or to knowledge of the 
cell temperature and laser power. For example, if the 
temperature were ten degrees lower than indicated, the 
limit for the 745 keV transition would become 8 > 11.1. 
With a change of temperature in the model, the best-fit 
value of r§° changes significantly to offset changes in the 
temperature-dependent values of the rubidium polariza- 
tion and the spin-exchange rate. The best fit (S = oo) is 
plotted for the 745 keV line in Fig. [4j 

In summary, radon has been collected and polarized in 
a glass cell, the temperature dependence of anisotropics 
has been measured, and new constraints have been set on 
the multipolc mixing ratios for 209 At transitions. These 
results are significant for an EDM measurement being 
planned at TRIUMF. The ISAC facility is expected to 
produce 223 Rn at a rate of 10 7 /s, more than 1000 times 
the rates for this work. At this rate, after two half-lives 
of implantation in a zirconium foil, approximately 500 
million radon atoms will have been collected. The pre- 
cession frequency of polarized radon in an electric field is 
proportional to the EDM, d. The expected uncertainty 
in d, when measured using the anisotropy signal, is 
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(5) 



° 2ET 2 y A 2 (1-5) 2 A' 
where E is the magnitude of the electric field, T 2 is the co 



TABLE I: The four main 209 Rn gamma ray lines (from the 
Table of Isotopes) are listed together with the corresponding 
absolute intensity I 7 , spin transition, and multipole mixing 
ratio (5). 



herence time of the polarization, A is the analyzing power 
of the measurement (i.e., the change in the signal due to 
a change in polarization), N is the total number of pho- 
tons detected, and B is the fraction of those photons due 
to background. Using T 2 = 30 s, based on the uncoated- 
cell relaxation time, A = 0.1 and B = 0.01, which are 
consistent with this work, and E = 5 kV/cm 0], a sensi- 
tivity of better than 10~ 26 e-cm is possible with N = 10 12 
gammas. Assuming a factor of 400 greater sensitivity for 
223 Rn, this would improve sensitivity to sources of CP 
violation by an order of magnitude or more compared to 
the limits set by the 199 Hg measurements. 
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